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Two-fluxon state in an annular Josephson junction in the presence of external magnetic field is 
studied analytically, numerically and experimentally. We obtain an analytical expression for the 
potential of interaction between the fluxons moving at arbitrary velocities (without the use of the 
"nonrelativistic" approximation). Treating the fluxons as quasi-particles, we then derive equations of 
motion for them. Direct simulations of the full extended sine-Gordon equation are in good agreement 
with results produced by the analytical model, in a relevant parameter region. Experimental data 
qualitatively agree with the numerical results. 
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I. INTRODUCTION 

A magnetic flux quantum (fluxon) in long Josephson 
junction (LJJ) is a well-known physical example of a sine- 
Gordon soliton. Ring-shaped (annular) LLJs serve as the 
ideal setting to study the fluxon dynamics, as it is not 
perturbed by boundary conditions at edges, which is the 
case for linear LLJs^. Due to the magnetic-flux quantiza- 
tion in a superconducting ring, the number of fluxons ini- 
tially trapped in an annular junction is conserved. An ef- 
fective tool, which makes it possible to create an effective 
spatially periodic potential for a fluxon trapped in the an- 
nular LJJ, is external dc magnetic field directed parallel 
to the ring's plane^-. If is the angular fluxon coordinate 
along the ring, the effective potential is U{6) ^ HcosO, 
where H is the strength of the magnetic fleld. The min- 
imum of the potential is located at the spot where the 
fluxon's magnetic moment is directed along the external 
fleld. The dynamics of a single fluxon in the spatially 
periodic potential has been an object of intensive theo- 
retical and experimental investigations^iiSiSiLSiSiifliii. 

The objective of the present work is to study, both 
theoretically and experimentally, dynamics of two flux- 
ons with equal polarities trapped in the annular LJJ. The 
system is schematically shown in Fig.^ This problem is 
distinguished by the interplay of the above-mentioned ef- 
fective periodic potential acting on each fluxon and direct 
interaction (repulsion) between them. 

The interaction between two fluxons in LJJs was flrst 
studied theoretically by Karpman et al. (see Ref. (3 and 
references therein). These authors found an analytic ex- 
pression for the interaction force between fluxons in the 
case of a small relative velocity, which corresponds to the 
"nonrelativistic" approximation when the relative veloc- 
ity is much smaller than the limit velocity in the LJJ (the 
Swihart velocity). However, in the situation relevant to 
experiment, the latter condition is not met, in the gen- 
eral case. In this paper, we aim to develop an analytical 
description of the interaction valid in the general ( "rela- 
tivistic" ) case. The theory will be based on an asymptotic 
method for weakly interacting solitons in non-integrable 
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FIG. 1: The schematic view of an annular Josephson junction 
with two trapped fluxons; the magnetic fleld H is applied in 
the plane of the junction. Josephson tunnel barrier is shown 
by thick black line; in grey are shown superconducting elec- 
trodes, which are extended in the junction plane in order to 
feed the bias current. 



simulations and presentation of experimental results. 

The theoretical model for the annular LJJ in the ex- 
ternal magnetic fleld was proposed in Refi^. It is based 
on the following extended (perturbed) sine-Gordon (sG) 
equation for the superconducting phase difference be- 
tween electrodes of the junction: 



iftt — sin if — aipt + "f + h sm{qx) . 
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The analysis will be followed by direct 



Here x is the coordinate along the ring, which is normal- 
ized to the Josephson penetration depth Aj, the time t is 
normalized to the inverse plasma frequency uJq^, a is a 
dissipation coefficient due to the quasiparticle tunnelling 
across the junction, and 7 is the bias current density 
(which is necessary in order to compensate losses due 
to dissipation), normalized to the critical current den- 
sity jc of the junction. As usually, the bias current is 
assumed to be uniformly distributed along the ring. Fur- 
ther, q = 2tt/L, where L is the normalized circumference 
of the junction, and h is the strength of the external mag- 
netic field H, normalized by a sample-specific geometric 
factori2i&. If N fluxons are trapped in the ring, Eq. ^ is 
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supplemented by the boundary condition 

ip{L + x,t) = ip{x, t) + 2ttN. 

The paper is organized as following. The single-fluxon 
dynamics in the annular LJJ is reviewed in Sec. Ill Al In 
Sec. IIIBI the derivation of an effective force of interac- 
tion between two fluxons, valid in the general (relativis- 
tic) case, is presented. Result of numerical calculations 



are displayed in Sec. IIIII In Sec. IIII Al we analytically 
consider a special case of an ostensible resonance in the 
two-fluxon system. In Sec. IIVI we present experimental 
results for two fluxons trapped in an annular LJJ. 



II. THEORY 

A. The basic model 

In our theoretical approach, we assume, as it is usually 
done, that the fluxons are well separated from each other, 
ICi ~ C2I S> 1, where ^1,2 are coordinates of their centers. 
In this case the two-soliton state may be represented by 
a linear combination of two single-soliton solutions: 



where 
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+ 27r(n - 1) (3) 



is the single-soliton solution of the unperturbed sG equa- 
tion, and ^„ is the velocity of the n-th soliton. The last 
term in Eq. |(2Jl is an arbitrary shift of the background 
phase (which is always possible if it is a multiple of 27r), 
which is chosen for convenience in what follows below. 

Before discussing the interaction of two fluxons, we 
briefly recall known results for a single fluxon trapped in 
an annular LJJ. This fluxon may be considered as a quasi- 
particle obeying the well-known equation of motionii^: 
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It is tantamount to the equation of motion for a rela- 
tivistic pendulum in a lossy medium under the action 
of constant torque. This equation has solutions of two 
types. The first type gives rise to \^{t) \ growing infinitely. 
It describes progressive motion (rotation) of the fiuxon 
around the ring, with a nonzero mean value of the veloc- 
ity ^. Solutions of the second type correspond to small 
oscillations of the fluxon around the minimum of the ef- 
fective potential with the frequency 
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the average velocity being zero. This state exists if I7I is 
below the critical value, 
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In the presence of dissipation (a ^ 0), the oscillations 
are damped, and in the stationary state the fiuxon is at 
rest. On the other hand, the progressive motion remains 
possible if the dissipation is not too strong. The fluxon 
moves with the average velocity, which, in the flrst ap- 
proximation, is given by the McLaughlin-Scott formula^. 
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Equation j?)) determines the current- voltage characteris- 
tics of the junction with a single trapped fluxon. 

In the case of two trapped fiuxons, we assume that 
they are quasiparticles interacting with a certain force 
(see below), and all the forces in Eq. act on each 
fluxon separately. In this case, three different dynamical 
regimes are expected: (i) oscillations of both fluxons ( "O- 
O" regime, i.e., a zero- voltage state) , (ii) rotation of both 
fluxons ("R-R" regime), (iii) rotation of one fiuxon and 
oscillations of the other one ( "R-0" regime) . Note that in 
R-0 regime oscillations take place even in the presence of 
dissipation, because the fluxon whose average velocity is 
zero is periodically excited by collisions with the rotating 
one. All the regimes are implied to be quasi-stationary 
states due to the stabilized effect of the dissipation. 



B. The interaction force 

In order to find the interaction force between two flux- 
ons, we use the center-of-mass reference frame (C-frame). 
In this frame, the two solitons move with velocities ±m, 
where 



1 - ^1^2 



^1-^2 

In the subsequent calculation, we set ^ < ^2- 

To calculate an effective potential of the interaction 
between fluxons, which will then produce the interaction 
force, we start with the Hamiltonian of the unperturbed 
sG equation for the infinitely long system: 



1 — cos ip 



(8) 



For the calculation of the Hamiltonian ||SJ|, we divide the 
space into two parts. The left (right) part occupies the 
space from —00 (+00) to the midpoint between the two 
fluxons, a = (^1 + C2)/2. We perform the actual cal- 
culation for the left part only, as for the right part the 
calculation is quite similar. 
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In the left part, we substitute the solution as the lin- 
ear combination where Lp2 is considered as a small 
perturbation, as the two fluxons are assumed to be well 
separated. Then, the Hamiltonian is written, in the first 
approximation, as 



left 



Hi + (5fficft,i 



where Hi = 8/VT — V? is the Hamiltonian of the unper- 
turbed sG soliton, and the interaction term is of the first 
order with respect to the weak field (^2 , 
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left.int 



dx ((pi^x^2,x + ^iASP2,t + ^2 sin ipi) (9) 



(recall the subscripts x and t stand for the correspond- 
ing partial derivatives). Substituting the expression valid 
for a moving soliton, ipn^t = —^'fn,x^ and integrating by 
parts, we obtain 



is the velocity of the center-of-mass of the two-fluxon set 
in L-frame. In the case of equal velocities, the potential 
(I12|l reduces to the well-known result of Karpman et ali^. 

The potential H12|l gives rise to two forces acting on 
each soliton, due to the ring geometry of the system, 
which should be added to the individual forces in Eq. 
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where AX = |^i — ■^2! , and 8 is for the effective mass 
of the fluxon in the present notation. Equations Q and 
(|14|l describe the dynamics of the two-fluxon system in 
an annular LJJ. These equations of motion we solved 
numerically by means of the fourth-order Runge-Kutta 
method. 



^-ffloft.int = (1 - '"^)'^l,x9'2|-oo + 
a 

+ / fix(/52 + <Pi,tt + sin(^i) . (10) 



The integral term in Eq. (|10|l is zero as the bracketed ex- 
pression is the sine-Gordon equation proper (this way to 
nullify the integral terms is known in the general analy- 
sis of the interaction between separated solitons^^). The 
contribution to the first term in Eq. H10() at the left limit, 
X = —00, is zero too because both (fi^x and Lp2 decay 
exponentially at infinity. In order to calculate the con- 
tribution from the upper limit, we use the asymptotic 
form of ipi,x {'^2) at large values of x (which is also a 
known point in the general analysis of the soliton-soliton 
interactionii): 
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After the substitution of the expressions Hll|l into Eq. 
(|10|l and calculation of the nonvanishing contribution to 
the first term from the right limit, a; = a, and then get- 
ting back from the C-frame to the laboratory reference 
frame (L-frame) , we arrive at an expression for the inter- 
action potential for two moving fluxons in the infinitely 
long junction. 
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where the contribution of the right half space is taken 
into account and 
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III. NUMERICAL CALCULATIONS 

In order to verify the theory presented above, we 
checked numerical solutions of the quasi-particle equa- 
tions of motion against direct simulations of the full equa- 
tion (P). Here we present results for fixed values a = 0.02 
and L — 20, while the bias current 7 was varied in steps 
of 0.002. 

The general behavior of the system can be described 
as the following. While 7 increases form zero, both flux- 
ons originally stay pinned in the effective potential in- 
duced by the magnetic field, so that the voltage across 
the junction is zero. At a critical value of the current 7c, 
the system switches to the R-0 regime, in which one of 
the fluxons rotates, while the other one oscillates due to 
periodic collisions with the moving fluxon. This state is 
stable up to another critical point, 7 < 73. On the other 
hand, decreasing the bias current leads to a transition to 
the regime with both fluxons pinned (zero voltage) at a 
different value, 7 = 7r. 

At 7 > 7s , the system operates in R-R regime with 
both fluxons rotating. Further increase of the bias cur- 
rent does not change the state of the system, up to a 
large value of the current, at which the junction switches 
to the "whirling" (alias resistive) state, with uniform ro- 
tation of the phase in all the system. When decreasing 
the bias current, the system switches first to the R-0 
regime, and then to the zero-voltage state. 

Typical current-voltage (I-V) characteristics, which 
display all these states and transitions, are presented in 
Fig. [5] Points shown by dots correspond to the numerical 
solution of the full equation , while the lines depict so- 
lutions of the quasi-particle model based on Eqs. and 
(|14|l . As is seen, the analytical quasi-particle model mak- 
ing use of the expression (|14(l for the interaction forces is 
in good agreement with direct simulations. 

The comparison of the critical values 7c of the bias 
current, obtained from the full simulations and from the 
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analytical model, is shown in Fig. O A small differ- 
ence between them can be attributed to the fact that, in 
the 0-0 regime, the actual distance between the pinned 
fluxons is small, hence the assumption of far separated 
fluxons does not apply in this case. Indeed the numerical 
simulations show that the distance between the fluxons 
in static case changes, depending on the magnetic field, 
in range of 0.8 — 1.5. 

For small currents, the I-V curve of the R-0 and R- 
R regimes, found from the direct simulations of Eq. 
feature additional small steps, which are due to resonant 
generation of radiation by the fluxons moving in the pe- 
riodic potential^. 

The comparison of the other critical value of the bias 
current, 7s (which corresponds to the first step of the 
I-V characteristics), again as found from the direct sim- 
ulations and analytical model, is shown, versus the mag- 
netic field, in Fig. 21 At small values of the magnetic 
field, these dependencies agree very well. However, for 
h > 0.45 the curve generated by the direct simulations 
goes down with the increase of the field. Actually, in this 
region the parameters 7 and h are too large to apply the 
perturbation theory. With the further increase of /i, 7^ 
decreases until it becomes equal to 7c. For a still stronger 
field, the system switches from the 0-0 regime directly 
to the whirling state. 

On the other hand, the curve produced by the analyti- 
cal approximation continues to go up with the field until 
h = 0.5. For fields larger than 0.5, the system resonantly 
switches from the R-0 regime to the R-R one, but at so 
large values of the field the quasi-particle model based on 
the perturbation theory becomes completely irrelevant. 
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FIG. 2: The current- voltage characteristics of the long an- 
nular junction with two trapped fluxons found from direct 
numerical simulation of Eq. Q (dots), and from the analyti- 
cal model based on Eq. (|1J and 1141 (solid line). In this case, 
the magnetic field is fixed to ft = 0.3. 
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FIG. 3: The critical current 7c versus the magnetic field h, ob- 
tained from direct numerical simulations of the full equation 
Q (dots) and from the quasi-particle model (line). 
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FIG. 4: The dependence of the maximum current of the first 
step in the I-V curve, 7^, on the magnetic field h, as found 
from direct simulations of the full equation Q (dots), and 
from the analytical model based on Eqs. Q ^-nd 1141 (line). 

A. A resonance condition 

An noteworthy feature of the two-fluxon dynamics in 
the R-0 regime is a possibility of a resonance between the 
natural frequency of oscillations of the trapped fluxon, 
and periodic excitation due to its collisions with the 
rotating one. The resonance condition is obtained by 
equating the frequency of small oscillations up, given by 
Eq. jSJl, and the rotation frequency Wr = Q < £,r > , where 
S,r is the velocity of the rotating fluxon, which can be ob- 
tained from Eq. This yields 

, /7r2\ I 64^^2 , ^ 

' = ^-^"^H^Jr^ ^^(-7.- + 4a)^ ' ^''^ 

where 7res is the value of the bias current corresponding 
to the resonance. 
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This dependence for L = 20 and a = 0.02 is plotted 
in Fig. 13 by dots. The line corresponds to the critical 
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FIG. 6: Current-voltage characteristics of the annular Joseph- 
son junction with two trapped fluxons. 



FIG. 5: Critical current of two fluxon state found from nu- 
merical solutions of Eq. Q (solid line) and the current of 
expected resonance (dots) 

current of the two fluxon state found from full numeri- 
cal simulations of Eq. for the same length. It may 
be expected that this resonance would result in a reso- 
nant switching from R-0 to R-R regime and drop of the 
critical current of the first step 7s vs. magnetic field h. 
However, at low magnetic fields the system switches to 
the R-0 branch at higher current than the current cor- 
responding to the resonance. At large magnetic field it 
switches directly to R-R branch, because the total per- 
turbation (field -I- bias current) is too strong. Here the 
perturbation approach is not any more applicable; the 
numerical solution of the analytical model reflects the 
resonance behavior by a drop on 7s (/i) dependence (solid 
line in Fig. With increase of the junction length L 
the the intersection point of the curves "fc{h) and ji-csih) 
moves upward in bias current, while its dependence on 
magnetic field is very weak. 



IV. EXPERIMENT 

Measurements of the I-V characteristics of the two- 
fluxon state were performed in long annular Nb-Al-AlOj;- 
Nb junctions. Due to the magnetic flux quantization in 
a superconducting ring, the number of initially trapped 
fluxons (solitons) is conserved. Trapping of a magnetic 
flux in the junction was achieved while cooling the sam- 
ple below the critical temperature T^^ = 9.2 K of nio- 
bium, in the presence of a small bias current passing 
through the junction. The number of trapped fluxons 
was determined from the highes voltage of the highest 
resonant branche on the current-voltage characteristics 
(Fig. ^ . Experiments were performed by applying the 
bias current / from top to the bottom electrode of the 
junction and measuring the dc voltage generated due to 
the motion (rotation, in terms of the theoretical consid- 



eration) of the trapped fluxons. The results presented 
below were obtained for a junction with the mean diam- 
eter 100 /zm and ring's width 3 /xm. The circumference 
(length of the annular junction) in the normalized units 
was L = 28.5, and the effective loss parameter was esti- 
mated as a = 0.03. The measurements were performed 
at 4.2 K. 

At zero magnetic field, depinning of a single-fluxon was 
observed, as a switching from the zero voltage state to 
the single-fluxon step of the I-V curve (in the state with 
one trapped fluxon), at the current Ic that was smaller 
by a factor w 65 than the critical current for the same 
junction, measured without trapped fluxons. This fact 
indicates a high degree of homogeneity of the junction 
(a strong local inhomogeneity would give rise to a much 
larger value of the fluxon-depinning critical current). 

We measured I-V curves of the state with two trapped 
fluxons for different strengths of the applied magnetic 
fleld. One of these curves is shown in Fig. The ex- 
perimental curves qualitatively agree with our analysis 
illustrated by Fig. |21 Two branches of the I-V curve are 
observed. The flrst branch, at around 65 /iV corresponds 
to the R-0 regime and the second one at about 130 /iV 
to R-R regime. 

Figure shows the measured value of the switching 
current for the transition from the R-0 regime to the R- 
R one, versus the external magnetic field. The curve is 
quite similar to the calculated dependence "fs{h) (dots in 
Fig- El I which was obtained above from direct integration 
of the full sine-Gordon model 



V. SUMMARY 

We have reported results of theoretical and experimen- 
tal studies of two-fluxon dynamics in an long annular 
Josephson junction in presence of the external magnetic 
field. An analytical expression for the interaction force 
between two fluxons moving at different arbitrary veloci- 



6 




0.0 0.5 1.0 1.5 2.0 2.5 

Magnetic field (a.u.) 

FIG. 7: The experimentally found critical current for the 
jump from the one-fluxon step of the I-V curve to the two- 
fiuxon one, as a function of the magnetic field. 

ties has been derived. Solution of the system of two cou- 
pled quasi-particle equations of motion for the fluxons, 
taking into account the interaction force, demonstrates 
good agreement with direct numerical simulations of the 
two-fluxon state in the full sine-Gordon model including 
all the perturbation factors. 

Three distinct dynamical regimes of the two-fluxon 
state have been thus identified. First, both fluxons may 
be pinned in the potential induced by the magnetic field. 
There is some discrepancy in the prediction of the critical 
current, which destroys this static regime, between the 



quasi-particle model and direct simulations, due to the 
fact that the separation between the two trapped fluxons 
is rather small, while the analytical model assumes them 
to be well separated. 

In the second regime, one of the fluxons rotates around 
the junction, while the other one oscillates in the poten- 
tial well. For this (R-0) regime, the maximum current is 
very accurately predicted by the quasi-particle model, if 
compared to direct numerical results. For this case, we 
have also investigated the possibility of the resonant ex- 
citation of the trapped fluxon by periodic collisions with 
the rotating one. 

In the third regime, both fluxons rotate. The corre- 
sponding I-V curves found by direct simulations demon- 
strate several additional small steps, which are due to 
resonant generation of small-amplitude plasma waves by 
fluxons moving in the periodic potential, which has been 
analyzed earliei»i for the single-fluxon case. 

The method developed in this work to calculate, in 
the general case, the effective interaction force acting be- 
tween two moving solitons, may find application to other 
soliton-bearing systems. 
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